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Abstract 

We show that curvature perturbations acquire a scale invariant spectrum for any constant 
equation of state, provided the fluid has a suitably time-dependent sound speed. In order 
for modes to exit the physical horizon, and in order to solve the usual problems of standard 
big bang cosmology, we argue that the only allowed possibilities are inflationary (albeit 
not necessarily slow-roll) expansion or ekpyrotic contraction. Non-Gaussianities offer many 
distinguish features. As usual with a small sound speed, non-Gaussianity can be relatively 
large, around current sensitivity levels. For DBI-like lagrangians, the amplitude is negative in 
the inflationary branch, and can be cither negative or positive in the ekpyrotic branch. Unlike 
the power spectrum, the three-point amplitude displays a large tilt that, in the expanding 
case, peaks on smallest scales. While the shape is predominantly of the equilateral type in 
the inflationary branch, as in DBI inflation, it is of the local form in the ekpyrotic branch. 
The tensor spectrum is also generically far from scale invariant. In the contracting case, for 
instance, tensors are strongly blue tilted, resulting in an unmeasurably small gravity wave 
amplitude on cosmic microwave background scales. 



1 Introduction 



The observation that the large scale structure originated from a nearly scale invariant primordial 
spectrum is generally construed as evidence that the universe underwent a phase of quasi-de Sitter 
expansion. Given the limited number of early universe parameters constrained by observations 
thus far, however, it is prudent to keep an open mind and seek the most general mechanism for 
producing density perturbations consistent with the data. It is well-known that generating a scale 
invariant spectrum of curvature perturbations with a single scalar field in 3+1 dimensions requires 
either inflation or a contracting dust-dominated phase [U On the other hand, scale invariant 
entropy perturbations can be generated in a larger variety of cosmological scenarios, provided 
more than one scalar field is present. This has been shown to occur, for example, in the slowly 
contracting background of the new ekpyrotic cosmology [31 HI O El d E] or in the kinetic dominated 
pre-big bang cosmology [HI [ID]- (See [TT] and [12] for recent reviews of the ekpyrotic/cyclic and 
pre-big bang scenarios, respectively.) Entropy perturbations can then be converted into curvature 
perturbations, for instance if the curvaton mechanism [131 El US] is at work. These earlier studies 
all have in common a constant sound speed Cg for scalar field fluctuations. 

In this paper we show that allowing for a time-dependent sound speed greatly broadens the 
range of possible cosmologies. The curvature perturbation on uniform density hypersurfaces C, 
can acquire a scale invariant spectrum for arbitrary background equation of state parameter 
e = —H/H"^ = 3(1 + w)/2, as long as the sound speed varies appropriately. 

Although our derivation applies to any perfect fluid, we take as our fiducial model a single 
scalar field with a Lagrangian 

£^ = v^P(X,0) (1.1) 

which is a general function of the kinetic term X = —^g'^'^d^(f)d^(f) and the field itself. In this 
theory the perturbations around a homogeneous background are still adiabatic. However, the 
presence of a non-standard kinetic term allows for a background dependent "speed of sound" of 
the fluctuation^ 

where p = 2XPx —P is the energy density of the field, and the quantities above are calculated on 
the unperturbed solution. While our basic mechanism for generating a scale invariant spectrum 
only relies on a varying sound speed of the adiabatic perturbations, we will refer more specifically 



to the single scalar field model (1.1) for calculating non-Gaussianities. 

In terms of the rate of change of the sound speed, = Cg/cgH, exact scale invariance requires 
e, to be constant and to satisfy a relation between them. More precisely, we find two branches 
of solutions: 



^Analogously, higher order curvature terms can induce a sound speed of tensor perturbations [16]. 
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• Case I: Expanding universe with decreasing sound speed, such that = — 2e. 

• Case II: Contracting universe with increasing sound speed, such that = |(3 — 2e). 

In the expanding case this mechanism was first proposed by |17j (see also P^ITU] ). whereas the 
contracting branch of solutions is new. Unlike [18], who considered a universe with superluminal 
sound speed initially, here we restrict our analysis to sub luminal sound speeds throughout. (The- 
ories with superluminal propagation have been argued not have a consistent UV-completion [20].) 

In the special case = 0, these recover de Sitter expansion and dust-like contraction, respec- 
tively. Of course, small deviations from these conditions and/or mild time-dependence for e,es lead 
to small departures from scale invariance. In order for perturbations to be stretched outside the 
Hubble horizon, the equation of state must satisfy e < 1 in Case I, corresponding to inflationary 
expansion, or e > 1 in Case II, corresponding to decelerated contraction. While the expand- 
ing branch requires inflation, it is worth emphasizing that our mechanism is not constrained by 
slow-roll conditions and therefore encompasses a much broader range of inflationary cosmologies. 

The above cosmologies, while degenerate at the level of the power spectrum, can be dis- 
tinguished through their non-Gaussian signatures. To begin with, the sign of the three-point 
amplitude is generically negative in Case I and can be either negative or positive in Case II. 
More importantly, unlike standard inflationary models, here the near scale invariance of the two- 
point function does not extend to the three-point amplitude. The non-Gaussian amplitude is 
proportional to evaluated at horizon crossing, and therefore exhibits strong running with 
scale [2T] [22] . In particular, non-Gaussianity peaks on the smallest scales in Case I since Cg de- 
creases in time, whereas the opposite is generally true in Case II. The sign of the tilt therefore 
allows to distinguish between the expanding and contracting branches. 

The shape of the three-point function in momentum space is also characteristically different in 
the two branches. In the expanding case, the three-point amplitude peaks for equilateral momen- 
tum triangles ^3], as usual for models with small sound speed such as Dirac-Born-Infeld (DBI) 
inflation [2U 1251 ES] ■ In the contracting case, however, the non-Gaussian signal features a strong 
peak for squeezed triangles, corresponding to the so-called local form p7]. This may come as a 
surprise since local non-Gaussianities usually require two scalar fields, as in the curvaton [H], mod- 
ulated reheating [28] , pre-big bang pump field [10] and New Ekpyrotic mechanisms [G] |29l |30l |3T] . 

In fact, the striking difference in shape between Cases I and II traces back to the evolution 
of C, at long wavelengths. In the expanding branch, the growing mode of C is a constant, and 
non-linearities have therefore a finite time to grow until modes go outside the horizon. Maxi- 
mal amplification is achieved when all wavevectors have magnitude comparable to the Hubble 
radius, leading to the predominantly equilateral shape. In the contracting branch, on the other 
hand, C keeps growing on super-Hubble scales, thereby allowing non-linearities to grow well after 
modes have exited the horizon. The non-Gaussian amplitude therefore peaks when one of the 
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wavenumbers is small, corresponding to the squeezed-triangle limit. 

The behavior of ( on super-Hubble scales is closely related to the issue of stability of the 
background. In the long- wavelength limit, ( represents a mildly scale-dependent perturbation in 
the scale factor [321 [33]. Thus all Case I cosmologies are stable, since a constant ( can be absorbed 
in a spatial diffeomorphism. Meanwhile, the super-Hubble growth of ( in Case II indicates that 
these backgrounds are unstable. Remarkably, such tachyonic instabilities appear to be a ubiquitous 
feature of contracting mechanisms: both New Ekpyrotic Cosmology and the general class of 
contracting two-field models of [SB] display a tachyonic instability in the entropy direction. This 
by no means represents a show-stopper — it was shown in |37], for instance, that a pre-ekpyrotic 
stabilizing phase can bring the field trajectory arbitrarily close to the desired classical solutions for 
a broad range of initial conditions. Moreover, as the present work underscores, the super-horizon 
growth of C leads to striking features in non-Gaussianity. 

Another important observable are primordial gravitational waves. Since tensor modes are 
oblivious to the time dependence of the sound speed and only feel the background cosmology, 
their spectrum is in general not scale invariant. In the expanding case, the gravity wave amplitude 
peaks on large scales, corresponding to a red tensor spectral tilt. In the contracting case, gravity 
waves are peaked on small scales, corresponding to a strong blue tilt. 

Various phenomenological considerations constrain the allowed equation of state parameter. 
We have already mentioned the inflationary condition e < 1 in Case I, in order for perturbations to 
exit the Hubble horizon. However, this range of e of course also ensures that (moderate amounts 
of) spatial curvature and anisotropy become increasingly irrelevant in time. If we impose similar 
standards in Case II, we obtain e > 3 (or w > 1), corresponding to the ekpyrotic range. Indeed, 
ekpyrotic contraction is equally powerful as inflation in solving the homogeneity, isotropy and flat- 
ness problems of standard big bang cosmology [21 EH]- Furthermore, the tilt of non-Gaussianities 
cannot be too strong, for otherwise perturbations risk exiting the perturbative regime within the 
observable range (Case I) or lead to an unacceptably large amplitude on cosmic microwave back- 
ground (CMB) scales (Case II). These non-Gaussian considerations lead to e < 0.5 in Case I. In 
Case II, on the other hand, non-Gaussianities are generically much larger and constrain the model 
in the region of the — e parameter space plotted in Fig. |2} Finally, since the tensor spectrum is 
red in the expanding case, we must ensure that the gravitational wave amplitude on large scales 
is consistent with CMB bounds. This requires the slighly tighter bound: e < 0.3 (Case I). 



While our analysis applies to general scalar field models of the form (1.1), we also provide 
explicit microphysical realizations of the scale invariant conditions within the context of DBI 
effective lagrangians. Here we piggyback on ^9], where it was shown that any scaling solution 
for the scale factor and sound speed can be obtained by suitably choosing the warp factor /(0) 
and scalar potential V{(j)). Remarkably, the desired condition = — 2e is satisfied (at least in 
the regime of small sound speed) for the AdS warp factor and m^0^ potential originally studied 



by [21] • Thus scale invariance obtains in this case without requiring the slow roll condition e ^ 1. 
To conclude this introduction, we offer an executive summary of our key results: 



A scale invariant spectrum of perturbations can be obtained for any background cosmology 
with constant equation of state, by suitably choosing the time dependence of the sound 
speed. We find both an expanding (Case I) and a contracting (Case II) branch of solutions. 
We discuss this basic mechanism in Sees. |2] and [H 

In the expanding case, various phenomenological bounds constrain the equation of state 
parameter to e < 0.3 and are discussed in Sec. |4} Our inflationary range is much broader 
than allowed by the usual slow-roll conditions. 

Scale invariance requires e,es to be constant and to satisfy a specific relation (see above). 
Deviations from these conditions lead to a non-zero spectral tilt, which is derived in Sec. |5| 

Gravitational waves are far from scale invariant. The tensor spectral index, derived in Sec. |6j 
is negative in Case I and positive in Case II, corresponding to red and blue tilts, respectively. 

Relatively large non-Gaussianities are generated in both cases (Sees. [Tjandjs]). In Case I 
they are of the DBI-equilateral type and generically negative; in Case II they are of the local 
type and generically positive. 

A (nearly) scale invariant two-point function does not necessarily imply a (nearly) scale 
invariant three point function. In fact, our mechanism fails to protect against running of 
non-Gaussianities. In Case I the three-point function has a strong blue tilt; in Case II it 



generally has a strong red tilt (see Sec. 8.2). 



2 Scale Invariance from Time-Dependent Sound Speed 

Our key observation is that in the presence of a time-dependent sound speed Cg the curvature per- 
turbation evolves according to an effective cosmological background different from the actual scale 
factor. Thus, a scale-invariant spectrum of density perturbations can arise for any background 
equation of state parameter e = —H/H"^, for a suitably- varying c^. In this Section, we show that 
scale invariance can also obtain in the contracting case as well. 

At quadratic order, the action for the curvature perturbation ( for general speed of sound 
models is given by [iO] 

(2.1) 
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S 



Ml, 



dC 

dr 



where r is conformal time, and z is defined as usual by 2; = ayltjcs. This holds whether the 
matter is a perfect fluid or a P(X, 0) scalar field. Since is a general function of time, it is 
convenient to instead work in terms of the "sound-horizon" time dy = c^dr: 

S = Ml,^ ^\&yq^\e-{^Cf \ , (2.2) 

where ' = d/dy, and 

- r- ^^2^ Q^ 
q = y/CsZ = —— . (2.3) 



The virtue of this time redefinition is manifest: the kinetic term takes the standard form, as if 
Cs = 1, while the dependence on Cg has been absorbed in a modified background described by q. 

In terms of the canonically- normalized scalar variable v = qC,, the equations of motion for the 
Fourier modes are given by 

(fc'-^)t;, = 0. (2.4) 

It is well-known that this results in a scale-invariant spectrum if and only if q" /q = 2/y'^. Indeed, 
for the adiabatic vacuum choice, the mode functions are given by 

^ » -iky 



resulting in a scale-invariant Pv{k) ~ fc^f^ in the limit k —>■ 0. Therefore, there are two possibilities: 

Q ^ {-yy^ or (2.6) 

where y runs from —00 to 0. In the hmit of constant c^, these of course reduce to de Sitter 
expansion and dust-dominated contraction PfllS], respectively. But, as we now show, more general 
backgrounds can also do the trick for suit ably- chosen c^. 

To proceed, let us assume for simplicity that the background has constant equation of state: 
e = const. Similarly, as usual we can define a rate of change parameter Cs for the sound speed, 

(2.7) 



He. 



also assumed to be constant. (In Sec. [sj we will relax these assumptions when calculating the 
spectral tilt.) From the definition of ?/, it is then straightforward to show that 

a{y) ~ i-y)^^ ; c,(y) ~ {-y)^^ . (2.8) 

Substituting these into ( 2.3[ ), then from (2.6) we deduce that scale-invariance occurs either for 



-2e (Case I) , (2.9) 
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corresponding to an expanding background with decreasing sound speed, or 



= ^(3-2e) (Case II) , (2.10) 

corresponding to a contracting universe with growing sound speed (for e > 3/2). Case I was pointed 
out in [T71 UHl [in] (see also [32] for a realization in DBI models of inflation |21]), whereas Case II 
is new. As a check, with = this recovers de Sitter expansion (e = 0) or dust-like contraction 
(e = 3/2). In slow- roll inflation, e is by definition small, and scale invariance forces to be small 
as well. What the above results show, however, is that much larger values of e are allowed as well, 
provided is chosen appropriately. 



3 Amplitude of Perturbations and Stability 



From the explicit solution for Vk in (2.5), we can read off the long- wavelength amplitude for 

Cfc = Vk/q: 

This can be further simpified by noting that ( |2.8[ ) allows to solve for y in terms H, a and Cg. 
Substituting the result into the above amphtude, the ( power spectrum is readily obtained: 

In the quasi-de Sitter and nearly constant Cg limit, e ^ 1 and es <^ 1, this agrees with earlier 
results PO] . 

Although the above expression for applies both to Cases I and II, the time-dependence 



of the right-hand side is different in each case. In Case I, H'^/cg is constant — see (2.8) — , 
and, following standard conventions, it therefore makes sense to evaluate at horizon-crossing 
(denoted by "bars"): 

P«-ii±^^ (3 3) 



where we have substituted (2.9). 



In Case II, on the other hand, H"^ /cg ~ l/y^, and the amphtude for C, keeps growing outside 
the horizon. Moreover, unlike Case I, here the final amplitude depends not only on the dynamics 
during the scale invariant phase, but also on the details of whatever phase precedes the bounce. 
For definiteness, we will assume that super-horizon modes stop growing after the scale invariant 
phase. This is achieved, for instance, if Cg saturates to a constant value and the collapse carries on 
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with an equation of state w > 1, or e > 3. In this case, the final amphtude of the power spectrum 
is evaluated at the end of the scale- invariant phase: 

(II) ^ (1 + e)^ Hl^^ 



where now we have substituted (2.10) 



Of course this difference in time-dependence is intimately tied to the stability of the back- 
ground. In the infinite- wavelength or A; — ^ limit, classical solutions for C, are interpreted as 
homogeneous perturbations to the background scale factor [32t [33]: 

ds^ = a\T){-dT^ + e^^dx^). (3.5) 



From C = ^/q and the explicit solution (2.5), we see that the growing mode is C ~ ^/qV- In Case I, 
which includes standard infiation, g ~ and ( —>■ const, on super-Hubble scales |3H|35]. Thus 
these backgrounds are stable. In Case II, however, q ^ y"^, and the growth of C as — *• indicates 
that the background is not an attractor. 

This instability appears to be a ubiquitous feature of all contracting mechanisms for generating 
perturbations. In the two-field ekpyrotic mechanism, for instance, the entropy direction must be 
tachyonic [Sj [TJ [37]. This was shown to be true for arbitrary two- derivative action of two scalar 
fields with general contracting backgrounds [36j. And here we find that these conclusions also 
apply to our single-field mechanism. 

Nevertheless, this instability should not be construed as a show-stopper. In the New Ekpyrotic 
scenario, for instance, it is circumvented by including a pre-ekpyrotic stabilizing phase that brings 
the field arbitrary close to the tachyonic ridge [37]. More importantly, as we will see later on, the 
growth of ( on super-Hubble scales enhances non-linearities and leads to interesting non-Gaussian 
signatures. 



4 Phenomenological Considerations 

In this Section we describe various constraints on the scale invariant phase, ranging from requiring 
that the modes, generated inside the sound horizon, cross the Hubble horizon, to solving the 
standard problems of big bang cosmology, to generating acceptable levels of non-Gaussianity. We 
will see that these conditions together require the background evolution to be infiationary (in the 
expanding case) or ekpyrotic (in the contracting case). 

4.1 Generating super-Hubble modes 

For general e, the amplification and freeze-out mechanism described above relies on a rapidly 
changing sound horizon. And indeed it is conceivable that the Hubble radius and sound horizons 
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are vastly different during this early-universe phase. But once the universe exits this phase and 
becomes radiation dominated, we must ensure that we are left with scale-invariant perturbations 
beyond the Hubble radius at reheating. In other words, during the generation of perturbations, 
modes must not only exit the sound horizon, but also the Hubble radius. 



Therefore the comoving Hubble radius, H ^/a, must be shrinking in time. From (2.8), this is 
given by 



-1 



(4.1) 

a 

and thus we require that the exponent be positive. This means e < 1 in Case I, corresponding to 
an inflationary universe, or e > 1 in Case II, describing a contracting universe with w > —1/3. 
While it is of course not surprising that the expanding branch must be inflationary, we should 
emphasize that this is not necessarily slow-roll inflation — here e can take on must larger values 
and still satisfy constraints on deviations from scale invariance. 

4.2 Duration of Scaling Behavior 

Although it is traditional to require scale invariant perturbations from the largest observable scales 
all the way down to the reheating scale, in reality observations only constrain near scale invariance 
over approximately 10 e-folds of modes, from CMB scales {k~^ ~ 10^ Mpc) down to galactic scales 
{k~^ ~ Mpc). Thus a conservative requirement on the duration of our scaling phase is 

Inf^^^M-L^ > 10. (4.2) 



aH/c 



sICMB 



As mentioned in Sec. |3| for instance, the scaling phase could be followed by a phase of nearly 
constant Cg, which would result in a different spectral tilt on sub-galactic scales. 

4.3 Homogeneity, Flatness and Isotropy 

We can also constrain e by requiring that the background solves the flatness, horizon and isotropy 
problems of standard big bang cosmology. Consider the general Friedmann equation 

2 K 8ttG I Cdust , Crad , Caniso \ , StcG 



where on the right-hand side we have allowed respectively for spatial curvature, dust, radiation 
and anisotropic stress contributions, as well as our general component with equation of state 
parameter e. 

In an expanding universe, the most dangerous term is spatial curvature — this is the flatness 
problem. Hence, our extra component must redshift more slowly than \jc? to be an attractor of 



the cosmological evolution. By inspection, this requires e < 1, or w < —1/3, which is again the 
condition for inflation. 

In a contracting background, on the other hand, the worrisome component is the anisotropy, 
which blueshifts as 1/a^. Thus a contracting universe typically becomes increasingly anisotropic, 
which is at the root of chaotic mixmaster behavior [HI |12]. In order to blueshift faster and 
win over anisotropy, our component must satisfy e > 3, or > 1 [21 |3H]. This is the ekpyrotic 
branch [3|, HH HH HHl HHl S?]. See JTl\ for a review of the ekpyrotic/cyclic scenario. Indeed, 
ekpyrotic contraction does equally well at solving the flatness, horizon and isotropy problems as 
inflationary expansion does. 

Two points are worth emphasizing: 

• The model with varying sound speed described here is the first example of a single-field 
ekpyrotic model resulting in an unambiguously scale invariant spectrum for the curvature 
perturbation. The earliest renditions of the ekpyrotic/cyclic scenarios also invoked a single 
scalar field, but ( was not scale invariant during the contracting phase, and one had to invoke 
higher-dimensional effects at the singularity [IHl HHl HHl |50] or a breakdown of the equation 
of motion for ( at the bounce [51]. In New Ekpyrotic Cosmology, this issue was resolved by 
using two ekpyrotic scalar fields to generate a scale-invariant entropy perturbation [31 IH El 

El Eli. 

• These earlier mechanisms, either relying on matching conditions at the singularity or with 
two scalar fields, all require e ^ 3 for the spectrum to be nearly scale invariant. Here the 
condition that e > 3 is much weaker, thereby extending the range of allowed models. 



While the range e > 3 is certainly desirable as it allows the ekpyrotic component to dominate 
over all other forms of energy density, this condition is by no means necessary from the strict 
point of view of generating perturbations. The actual lower bound comes from generating super- 
Hubble modes which, as we have seen in Sec. 4J^, only requires e > 1. For instance, the instability 
in Case II presumably requires some earlier phase to bring the field trajectory sufficiently close 
to the desired classical trajectory; this pre-scaling phase could also drive the universe towards 
homogeneity and isotropy using the ekpyrotic smoothing mechanism if e > 3. 



4.4 Constraints from Gravity Waves and non-Gaussianities 

Further constraints to our model come from the tilt of the gravitational wave spectrum and non- 
Gaussianity, which will be discussed in more detail in Sees. |6] and 8.2[ respectively. Interestingly, 
the cancelation between e and that results in a scale invariant two-point function does not 
protect the tensor spectrum nor the three-point function from running. 
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The inflationary branch (Case I) is characterized by a decreasing sound speed, which in turn 
leads to larger non-Gaussianities on small scales. In order to maintain perturbative control over 
the entire range of modes, we must require e < 0.5 (or w < —2/3). Meanwhile, the tensor spectrum 
is red and therefore peaks on large scales. To satisfy current bounds on the gravitational wave 
contribution to the CMB temperature anisotropy, we must impose e < 0.3. 

In the ekpyrotic branch (Case II), the growing sound speed results in non-Gaussianities that 
grow in amplitude as we probe larger scales; but clearly the sound speed should not exceed the 
speed of light during the evolution. Moreover, the growing of the modes outside the horizon 



generate a new plethora of terms that have to be included in the three-point function (see (7.39)) 



and that can very easily produce unacceptably large level of non-Gaussianity on CMB scales. 



However, we can still find a safe region in the (cc,end, e) parameter space. With /x = 1 — see (7.33) 
for its definition — , for instance, the allowed region lies between Cg > 0.7 and 2.8 < e < 3.15, as 
shown in Fig. |2} Meanwhile, the tensor spectrum has a strong blue tilt and therefore does not 
yield any constraint on the equation of state parameter. 

To summarize, for our two branches of solutions, corresponding respectively to expansion and 
contraction, we have the allowed respective range for the equation of state parameter: 

< e < 0.3 (Case I) 
2.8 < e < 3.15 (Case II). (4.4) 

Intriguingly, the bounds on the equation of state are comparably constraining in each case: roughly 
up to 20% allowed deviations from exact de Sitter (—1 < w < —0.8) or kinetic domination 
(0.9 < w < 1.1), respectively. Note that in Case II the sound speed must be suitably adjusted 
as a function of e in order to give acceptable non-Gaussianity, as shown in Fig. [2] Each branch 
addresses the standard problems of big bang cosmology and, for suitably varying speed of sound, 
yields a scale invariant spectrum of density perturbations while being consistent with present 
limits on gravitational waves and non-Gaussianity. 



5 Scalar Spectral Index 



In this Section we consider small departures from exact scale invariance, which occur whenever 



deviates from the conditions (2.9) or (2.10), and/or when e, are time-dependent. To parametrize 



this time dependence, we introduce 



1 dine 
~H dt 



1 dines 
H &t 



(5.1) 



For simplicity, we treat rj and rj^ as small and constant. Nothing else is assumed about e and e^ 
for the moment. 



10 



In this approximation, our time redefinition, dy = c^dr, can be integrated explicitly: 

er] + esiis 



y 



(e + e, - l)aH 



1 + 



[e + ts 



(5.2) 



And this in turn can be used to obtain an explicit expression for q" /q in the mode equation (2.4). 
For instance, we have 



ding l + {rj-es)/2 ( erj + esVs 



(5.3) 



dlny e + e, -1 ' (e + - 1)1 

Proceeding along these lines, after some straightforward algebra the mode equation takes the form 



where 



1 



1 



4 {e + ts-lf 

1 + - e,)/2 
e + e, - 1 



^ V-(^s \ esVs {1 - es/2){A- 3es - 2e){eT] + esT]s) 
' 2 + (e + es-iy 



(5.4) 



(5.5) 



As usual we can turn (5.4) into a Bessel equation for Vk/\/—ky, which, assuming adiabatic 
initial conditions, has a Hankel function solution, Hl^\—ky). In the long-wavelength limit, the 
corresponding spectral index is then 

-l = -2i^ + 3. (5.6) 
Thus the spectrum is nearly scale invariant for z/ ^ 3/2. In particular, setting rj = ris = 0, one can 



quickly check that exact scale invariance {u = 3/2) requires that either (2.9) or (2.10) be satisfied. 
We next derive explicit expressions for for small departures from these conditions. 

• Case I: ^ — 2e 

It is convenient to characterize deviations from = — 2e through a small parameter 6i: 

e, = -2e + (l + e)5i. (5.7) 
By expanding p to linear order in 6i, r], and rjs, the scalar spectral index is thus given by 



K - i)i 



-5i 



1 



;i + e)2 
11 



. 11 \ 14 ■ 



(5.^ 



As a check, note that in the hmit of quasi-de Sitter background and nearly constant c^, for 
which e, -C 1, we have (5i ~ + 2e, and thus 



{Us - 1) 



slow— roll 



'2e- es-Tj . 



(5.9) 



This agrees with earher treatments [40J. 

More generally, however, larger values of e and e<j are possible without spoiling approximate 
scale invariance, provided 5i remains small. Thus the 5i parametrization offers a more general 
characterization of deviations from = 1, which, in the special case e, eg -C 1, reproduces 
standard slow-roll results. 

Instead of 77^, we can alternatively introduce a parameter 5j obtained by taking the derivative 

d5i 



of (5.7): 



5' 

^ d log a 

Expressed in terms of this parameter, we find 

(n, - 1)t = -di - 7:;— — 



n 



1 + e 



31 + e 



(5.10) 



(5.11) 



Case II: e, ^ §(3 - 2e) 

Similarly, in the contracting branch we can introduce Su as 



e, = ^(3-2e)-^(l + e)5„ 



(5.12) 



Substituting into (5.5) and, as before, keeping terms linear in Su, r], and r]s, we obtain 

5 



K - l)ii = -5 5ii - ^^^^y 

Or, equivalently, in terms of 5jj = d(5n/dloga 

5 

{us - 1)11 = -5 5u 



'l + e] 



43 \ / 2 

1 + ye ) ^ + 70 n--e)r/. 



^ 97 \ 35,^ 



(5.13) 



(5.14) 



In slow-roll inflation e and t] are usually considered as independent and comparably small 
parameters. Similarly, here we treat 6, t] and 77^ (or 6, rj and 6') on the same footing because 
that all appear at linear order in the expression of the observationally small quantity n<j — 1. Of 
course, more fine tuned cancellations and richer hierarchies between the parameters can also be 
considered. 
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6 Gravitational Wave Spectrum 



Unlike scalar perturbations, tensor modes are oblivious to a time-dependent sound speed and 
instead feel only the cosmological background. The mode functions hk{T) therefore satisfy 

As usual, this can be cast as a Bessel equation through the field redefinition: Uk = ahk/ ^/—kr. 
Assuming adiabatic initial conditions, the solutions are given by 

hkiT) = ^^H^^l i-kr). (6.2) 

• Case I: ~ — 2e 

In the inflationary branch, the equation of state parameter falls in the range < e < 1. The 
tensor power spectrum, Ph = /c'^l/i/cp/vr^, therefore takes the form 



H^„^j ( k 



2e 
" l-£ 



corresponding to a red tensor tilt: 



n?' = ^^ = -^. (6.4) 

d m K 1 — e 



The tensor-to-scalar ratio follows by comparing with (3.3). Noting that the sound speed 



evaluated at horizon crossing satisfies ~ /c we obtain 

r« ~ ec. (^^J . (6.5) 

And since this amplitude peaks on the largest scales, we have to ensure that this satisfies 
current bounds on the gravitational wave contribution to the CMB: tcmb ^ 0.3. We will 
see in Sec. [sjthat the non-Gaussian amplitude is proportional to c^^, and thus the current 
constraint from non-Gaussianity gives Cs|cmb ^0.1. Moreover, as discussed in Sec 



4.2 



we 



have i^endZ-f^lcMB ^ 10^. Assuming the minimum required number of e-folds, we obtain the 
conservative bound 

e<0.3. (6.6) 



This confirms the bound quoted in (4.4) 
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• Case II: e, ^ §(3 - 2e) 

In the ekpyrotic branch, the gravitational wave power spectrum is 



,(II) -^end 



2 / h \3-lf5fl 



^ Ml, V^end 

For e > 3 the tensor spectrum has a strong blue tilt: 



(6.7) 



4") = -^>2. (6.8) 

In particular, the gravitational wave amplitude is negligible on CMB scales in this 
in the pre-big bang [9l [52] and ekpyrotic scenarios [ISj [53]. For e < 3 we have 

= (6.9) 

The gravitational wave spectrum is flat for e = 3/2, i.e., in the case of dust-like contrac- 
tion pLj. 



7 Calculating Non-Gaussianities 



We now calculate the three-point function for the following general single scalar field model [261155] : 



S = d XV— ^ 



2 + 



(7.1) 



where the pressure P is a general function of the scalar field and the kinetic term X = 
— ^g'^'^d^(f)d^(f). To simplify the expressions below, we set Mpi = 1 throughout the calculation. 
The energy density reads 

p = 2XPx-P, (7.2) 

while the speed of sound is given by 

Besides the parameters e, e^, r], rjs introduced earlier, the calculation of the three-point function 
also necessitates defining two further parameters derived from P{X, 0) [55[ 126] 

E = XP,x + 2X''P,xx = ^ , (7.4) 

A = X'Pxx + \x^P,xxx . (7.5) 
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Consider for instance the DBI effective action [21] describing the motion of a 3-brane in a 
throat of some higher-dimensional compact space: 



-Pdbi(-^, 



(7.6) 



where /(</>) encodes the warp factor of the throat (see e.g. [Sj for non-standard throat geometries). 
In this case we have 



c^ = l-2/(0)X; 



X 



;i - 2/(0)x)3/2 



A 



XV(0) 



;i - 2/(0)x)5/2 • 



(7.7) 



Although our calculation of non-Gaussiniaties below applies to general P{X, 0), we view the DBI 
action as a fiducial model for our analysis. In particular, in Appendix |A] we follow ||39| and provide 



explicit examples of DBI Lagrangians that yield scaling solutions of the form (2.8) 



In the present case, the cubic effective action has already been derived in [261 [55] . The result 
is vahd outside of the slow-roll approximation and for any time-dependent sound speed: 



(3) 



dtd^x 



S I 



1 



2A 



+ 



3 + 3c^)CC^ 



ae , 



ci)adcf-2a-adO{dx) 



+ ;r^^ I ^ I C'C + Ya^dO{dx)d\ + {-^{d\){dxf + 2/(0 ^ 



e 

4a' 



5L( 



(7.1 



where dots denote derivatives with respect to proper time t, d is a spatial derivative, and x is 
defined as 



-C- 



(7.9) 



5L, 



Meanwhile, in the last term -^|i denotes the variation of the quadratic action with respect to 
the perturbation C,: 



6L 



(2) 



(7.10) 



/(C) 



+ 



+ —a + — 



m){dx) - d-\d,d,{d,cd,x))] 



(7.11) 



where (9 ^ is the inverse Laplacian. In order to calculate the three-point function in the interaction 



picture we can derive the Hamiltonian H^^x, in the usual way. Except for the last term in (7.8) 
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we have ifjnt = ~L{3), up to interactions that are higher order in the number of fields. The last 
term requires more attention because it contains a second time derivative in the field which makes 
problematic the definition of the conjugate momentum. However, since f^|i is proportional to the 
linearized equations of motion, it can be absorbed by a field redefinition 

C-Cn + /(Cn). (7.12) 



Only the first two terms in (7.11) contribute to the three-point correlation function at long 
wavelengths since all other terms involve gradients that go to zero outside the horizon. In Case I, 
the second term is also negligible since ( — > const, on super Hubble scales. Thus, after the field 
redefinition C Cn + ^^Cn' three-point function can be expressed in this 

(C(X1)C(X2)C(X3))(I) = (Cn(xi)C„(x2)Cn(x3)) 

+ ^ ((C„(xi)Cn(x2))(Cn(xi)Cn(x3)) + Cyclic) + 0{C) ■ (7.13) 

In Case II, on the other hand, the second term can be simplified by noting that C ~ l/y'^ on 



super-Hubble scales. To leading order in ri,ris, it follows from (5.2) and the condition for scale 



invariance (2.10) that C„ = — (3/5)if(l + e)^„. Thus the desired field redefinition is C — ^ Cn + 
i^Cn ~ ^^5^'^'* Cn; tcrms of which the three-point amplitude takes the form 

(C(X1)C(X2)C(X3))(„) = (C„(Xi)Cn(x2)Cn(x3)) 

+ - ^^^) ((Cn(xi)Cn(x2))(C„(xi)C„(x3)) + Cyclic) 

+ OiC). (7.14) 



7.1 The three-point function: some useful formulae 

We now proceed to calculate the three-point function using standard methods [33j . At first order 
in perturbation theory and in the interaction picture we have 

(C(t,ki)c(t,k2)c(t,k3)) = -z f dt'{[atM)atM)atM),Hut')]) , (7.15) 

J to 

where to is some sufficiently early time. It proves convenient to work with the sound horizon time 



variable y defined in (5.2) [dy = Csdt/a) and expand the field ( in creators and annihilators in the 
usual way: 

Ciy, k) = Ukiy)aik) + <(|/)at(-k). (7.16) 

Note that with our Fourier conventions commutation relations read [a(k), a'''(k')] = (27r)^(5^(k— k'). 
The mode functions Uk are easily computed in terms of the canonically normalized mode functions 

16 



V (2.5). Up to an irrelevant phase, we have 



Uk{y) 



C \l/2 



Hil 



'1 + iky)e 



-iky 



(7.17) 



The contribution to the three-point function of each term in (7.8) can be worked out once we 



substitute (7.16) and (7.17) in (7.15). After switching to the integration variable y and making 



some simplifications, we encounter integrals of the form 



C 



Vend 



dy 



Vend 



(7.18) 



Only the imaginary part of this integral is relevant for our calculation. In the above expression, 
n is an integer, y^nd < is the value of y at the end of the inflationary or ekpyrotic phase, and 
K = ki + k2 + k:i. Hence, for all modes of interest, K\yend\ is a small quantity. As discussed in [33] . 
the above choice of integration contour picks up the appropriate interacting "in-in" vacuum at 
\y\ ^ oo and takes care of the oscillating behavior of the exponential. 

For simplicity, we consider only strictly scale invariant solutions and therefore neglect small 
corrections of order 5, rj and rjs in our calculations. The exponent 7 is thus a constant, parame- 



terizing the behavior of either HjcJ'^ or H/cJ^ ■, depending on the term being considered. More 
specifically, in Case I (expanding branch) we have 



^5/2 



H H, 



end 



.1/2 



-1/2 
end 



(7 = 0) 



H H, 



end 



where 



In Case II, instead. 



5/2 5/2 , ^ 



4e 



(7 



a 



a 



1 + e 



H H, 



end 



.1/2 



-1/2 
'^s end 



Vend 



-3 



(7 = -3) , 



H H, 



end 



^5/2 



.5/2 
''s end 



Vend 



where 



p = 5- 



(7.19) 



(7.20) 



(7.21) 



(7.22) 



Let us describe the integrals in more detail for Case I. For 7 + n > —2 the imaginary part 
of (7.18) is convergent as y^nd 0. In this case we can approximately extend the upper limit of 
integration to 0, which amounts to neglecting terms of higher order in fc|?/end|- We thus obtain 



ImC 



■(irlyendl)-^ COS ^r(l + 7 + n)K-^-' 
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(7.23) 



When 7 + n < —2, however, the integral (7.18) is divergent for i/cnd 0. To leading order in 
k\yend\ we find 



ImC = Im- 



7 + n + 



(7.24) 



For example, 



ImC = —— m = 

i^(7 + 2) ^ ^ 



ImC 



'i^2(^ + 2) 



(n = 1) , 



(7.25) 



and so forth. For the parameter range relevant to Case II the integrals are calculated in Appendix 

m 



7.2 Example 1: the (('^ contribution 

Before writing our full set of results, it is worth outlining as an example the calculation of the 



3-point contribution from the (("^ term in (7.8). By applying the commutation relations this reads: 



(C(ki)C(k2)C(k3)) 



end I 



X 



Vend 



dy 



Cs a e 



e-3 + 3ci)ul{y) 



dy dy 



+ perm. + c.c. (7.26) 



Substituting (7.17), we have 



(C(ki)C(k2)C(k3)) 



t{27ry6'{ki + ks + kg 



^id(l-e-e.)3 1 



43.3/2^3/2 111.3 
^ ^ '-send ^ ■? 



"J/cnd 

X / dy 

oo+ie 



— (e-3 + 3cJ ^3/2^3/2 



[I - ikiy)klkly'^e'^y + perm. + c.c. (7.27) 



Since we are considering a strictly scale-invariant solution, e and are constant and can be taken 
outside of the integral. Now, we can use (5.2) to substitute for a^y"^ and find 



(C(ki)C(k2)C(k3))c^2 = i{2ixf5\k, + k, + W 



X / d?/^(e-3 + 3c^)(l-iA;i?/)A;^fc^e^^^ + perm. + C.C. (7.28) 



oo+ie C, 



5/2 



In the expanding case (Case I), 



-2e, and the background-dependent functions in the 



integrand display the two types of time dependence listed in (7.19). Since both of these give 
rise to convergent integrals, we can apply (7.23). Moreover, we combine the prefactor {K\ycnd\)~'^ 



with the appropriate background functions, Hcnd/cJ^nd -^end/c^end' express the corresponding 
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J>/2 



quantities at horizon crossing. In other words, H^ndC 

-5/2 



"1/2 
s end 



5/2/ 



Hcs where bar quantities are evaluated at horizon crossing. We finally obtain (Case I) 



(C(ki)C(k2)C(k3)) 



CC^ (I) 



X 



{2Txf5\\ii + k2 + k3 



3) cos —T[l + a) 



'^2 '^3 



16e2cf Ylj^ K 



1 + (1 



a 



K 



3c 



7,2 



K 



sym. (7.29) 



Note that the combination cos ^r(l + a), which appears ubiquitously in Case I (see below), is 
regular in the limit a ^ — 1: the constraint « ^ — 1 (e ^ 0.3) is uniquely dictated by phenomenol- 
ogy- 

In the contracting case (Case II), = 2(3 — 2e)/5, and the background-dependent functions 



in the integrand follow the power-law behaviors given in (7.21). Moreover, since the amplitudes 



of the perturbations keep growing outside the horizon, it is convenient to refer to the background 
quantities at the end of the ekpyrosis phase as those directly relate to the power spectrum. This 
growth of the amplitude complicates the calculation somewhat compared to the expanding case. 



Specifically, the time-derivative of the mode function (7.17) in Case II reads 



dy 



k'^y (1 - iky) 

y 



(7.30) 



Note that only the first term in square brackets appears in Case I. The presence of the extra terms 



results in longer expressions for the amplitudes in Case II. The final answer is given in Sec. 7.4 



7.3 Example 2: the contribution 

As another example of our method, consider the term, for which the three-point function reads 



(C(ki)C(k2)C(k3)) 



-i{2nf5^iki + k2 + k3)Uk^{ycnd)Uk2{ycnd)Uks{ycnd) 

'a2 i/sl 



C3 

X / dy— — [S(l-c^ ) + 2AJ — hperm. + C.C. (7.31) 



dy dy 



dy 



Note that the parameter A can be written as [55] 



6 V 



where we have introduced 



ee.. 



X 



3e 



X 
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(7.32) 



(7.33) 



with tx the "kinetic part" of e: 

e = e</, + ex 



dH 



X dH 



ex 



(7.34) 



For simphcity we will henceforth assume that ex is constant. Note that for the DBI action ( 7.6[ ) 

/r = 1 - (7.35) 

which therefore lies between and 1. It follows that the contribution vanishes identically for 
the DBI case [26] . 

Explicitly, in the expanding case, we obtain 



(C(ki)C(k2)C(k3))^3 (I) = -(27r)V(ki+k2 + k 



an 



2c' + (/x - l)cos— r(3 + a) . (7.36) 



Again, for Case II, we refer to Sec. 7.4 



7.4 Combining all terms 

The remaining contributions are calculated in a similar way to the above examples. After com- 
bining all the terms, it is convenient to express the full three-point function by factoring out 
appropriate powers of the power spectrum and defining an amplitude A as 



(C(ki)C(k2)C(k3)) = {2ny6%k, + k2 + ks)P^'—^A., 



(7.37) 



where is given in (3.3) and (3.4) for Cases I and II, respectively. Collecting all terms, we obtain: 



Case I: 



A 



(I) 



1 + e 
2c2 



an ^ 



^ci + {fx - 1) cos — r(3 + a 



A 



(I) 

4c2 



6cl + (e - 3) cos ^(2 + a)r(l + «) I ^ ^ k^k^+ 



K 



i<j 



3c1 + (e-3)cos^r(2 + a)l-l5^fc,'fc| 
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A 



(I) 



(1 + 5e) Q!7r„, , 
^ ' cos — r(l + a) 



(I) 



cos — r(l + a) 



2 ^ -J - ^2 



"^^^ 16c2 



cos— r(l + Q;)(2 + Q;/2) 



.2,, „„J ('a2A\/Q,A2 



X2 



(7.38) 



where A^^ accounts for the d(dxd x (9 0(dx) terms. 



• Case II: 



(11) 



A\' 



7r/3 
2 



(l + e)(/x-l) [ cos ^ 

"\(i^l2/end|)^ 



^'-s end 



r(/3 + 3)M|M + 3r(/? + 2) 



+ 



9K' 



(/3-2)(/3-3) 



cosf ^^^^^M±3^3 



(3-1 



- 3(1 + e)K'^ 



A 



(11) 



(e-3) I cosf 

4cLnd l(i^lyendl)^ 



r(/3 + i) 



2 + /3 



+ 3r(/3) ( 2 5^ hk] + 9A;iA;2A;3 + 



i<j 



(/3-2)(/3-3) 



i<j 



i<j 



3 (3^-5(3-3 
(/?-2)(/3-3) 
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1 +/? 



A 



A 



(11) 

2 



+ 



[ 2K ^ hkj + SX^ + ^ hk] 



(II) _ 13e-7 cos^ 



r(/5 + l) J^A:,^ 



K 



^'^s end 



cosf r(/j + i) 



^ - 2{f3 + 3)kik2k3 



( COS f r(/5) 



2Kkik2k3 



+ Sir 



27 



cos 



nl3 



r(/3) , 



(i^bend|)'^/3-2^/?(/5-2) 3 



(3-3 



32C^ end 



COS fr(p + 1) 

(i^|?/end|)^ 



i2k,k2k, -^y: ^^^I + ]^ (E ^.^^J - 2 E f^f^ + E ^^^1) 

i¥=j i¥=j ij^j 



i+3 

/cosf r(/3) _ l' 



-5Akik2k3 - 15 kik^ ~ A;| + 



kik2k3 



3E* 



i^j i^j i<j 

cosfr(/3-i) 1 



E*?*? + E*?*J-2E*?*? 

. ^T^i *7^i i¥=j 



+ 3K 



- E^^ 



{K\yen,\r 
1 



6 ^ /c| - 18 ^ - 



i<j 



kik2kz 



"'1"'2"'3 



2E^^^-4E^^^+3E^4-3E^^^ 



i<j 

22 



+ 9 



cos 



7r/3 



(11) 

redef 



iK\yend\r/3 
3(l + e 



+ 



PiP - 2) 



f3-3 



6 - 




(7.39) 



where ^redef arises from the field redefinition (7.14). 

Note that the expressions in Case I and II take nearly identical form in the limit of small 
sound speed, under the replacement a ^ (3. The only significant difference is the overall factor 
of c^^ in Case I compared to c^^ ^^^{\ycnA\KY in Case II, which stems from the power spectrum 
being normalized at horizon-crossing and at the end of the contracting phase in Cases I and II, 
respectively. 



8 Non-Gaussian Signatures 

In this Section, we discuss three observable features of the 3-point function, each of which can 
distinguish the expanding branch from the contracting branch: the amplitude, parametrized by 
the /nl parameter; the spectral dependence, characterized by the tilt ^ng; and the momentum 
shape of the 3-point amplitude. 

Because of the reduced sound speed the 3-point amplitude is relatively large and potentially 
measurable by current or near-future CMB experiments. Interestingly, the /nl parameter is 
generically negative in the expanding case (including DEI infiation) and positive in the contracting 
case. In particular, this holds true for DBI lagrangians. Nevertheless, it is certainly possible to 
construct models where this conclusion is altered. 

In standard slow-roll infiation, including DBI, the quasi-de Sitter nature of the background 
(e, ?7 <^ 1) and the near constancy of the sound speed (cs ^ 1) imply both a nearly scale-invariant 
power spectrum as well as nearly scale invariant non-Gaussianities. In the generalized models 
studied here, however, the scale invariance of the 2-point function is maintained by balancing 
a relatively large e against a correspondingly large e^. This cancellation fails at the three-point 
level, leading to a strong scale dependence for non-Gaussianities. Here again the expanding and 
contracting solutions generically have a different signal: the inflationary branch (Case I) leads to 
a blue tilt for the 3-point amplitude, whereas the ekpyrotic branch (Case II) has a red tilt. 

But the most striking difference comes from the shape of the non-Gaussianities. In the in- 
flationary branch, the shape is predominantly of the equilateral type [23], which is characteristic 
of inflationary models with low sound speed, as in DBI [25]. In the ekpyrotic branch, however, 
remarkably the shape is predominantly local. 

This stark contrast in shape arises from the different behavior of C, on super-horizon scales. 
Because C, tends to a constant outside the horizon in Case I, non-linearities can only grow for 
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a finite period until modes exit the horizon. Thus non-Gaussianities peak when all wavevectors 
have magnitude comparable to the Hubble radius, and therefore comparable to one another — 
hence the equilateral shape. In Case II, ( keeps growing outside the horizon, thereby allowing 
non-linearities to grow until the end of the scaling phase. The non-Gaussian amplitude therefore 
peaks when one of the wavenumbers is small, corresponding to the squeezed-triangle limit. 

There is also an interesting subdominant contribution to the non-Gaussianities. While this is 
suppressed by slow-roll parameters in standard DBI inflation, here it accounts for a more significant 
fraction of the amplitude since e can be relatively large in Case I. In the inflationary branch, we 
find that the subdominant contribution is of the squashed shape; in the ekpyrotic branch, it is of 
the equilateral shape. 

8.1 Amplitude of Non-Gaussianities 

The 3-point amplitude is characterized as usual by the /nl parameter p7] : 

C = Q{x) + lhLCl (8.1) 

Here we follow the WMAP sign convention, where positive /nl physically corresponds to negative- 
skewness for the temperature fluctuations. The resulting 3-point function in momentum space is 
given by 

(C(kOC(k2)C(k3)) = {2ny6'{k, + k, + ^s)P'^h^^3 ^ (8-2) 

which peaks in the squeezed limit, e.g., ki <C k2, k^. 

Although the momentum dependence of our amplitudes is manifestly different, in particular 
peaking in the equilateral limit in Case I, it is conventional to define /nl by matching amplitudes 
at ^1 = ^2 = ^3 = K/?): 

/nl = 30 ^^-=^3^=^^ . (8.3) 

In the expanding case, by evaluating the right-hand side for the amplitudes listed in Sec. |7.4[ 
we obtain 

.(I) ^ 35 40 

108 243(a + 4) 

cos — r(l + a) { —{a + l)(a + 2) {fx - 1) + -(1 + a) + — ^ . (8.4) 



9c2(a + 4) 2 ' '[27' ' '6' '32 

To gain some intuition on the parametric dependence, a useful fitting formula for this expression 
is 

fSl ~ 0.27 - 0.164 _ (0.12 + 0.04/x) c^^l + a) , (8.5) 
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where we recall from (7.20) that —2 <a<OforO<e<l. Note that /nl is fairly insensitive to 
fx, for reasonable values of this parameter, including DBl. Figure [T] compares the fitting formula 
with the exact expression as a function of a and for various values of fx- The agreement is 
typically within 10 % over the phenomenologically allowed range of e < 0.3 (or a > — 1). 




Figure 1: Comparison between the exact formula (8.4) and its linear fitting (8.5) in the expanding 



case, for the parameter region of interest: — 1 < a < 0. Values of f^\ are evaluated as a function of 
a for different c, and fx'- from top to bottom, (c„ fx) = (0.25, 1); (0.15, 0.1); (0.12, 0.9); (0.1, 0.1). 



In the contracting case, on the other hand, the expression for /nl is much more involved, as 



can be seen from (7.39). Non-gaussianities are quite large for generic choice of parameters in this 
case. However acceptable values of /nl can be obtained in some region of the parameter space. 
Figure |2] shows the region in the (c^end^e) parameter space corresponding to —100 < /nl < 100, 
with /x = 1. 



8.2 Spectral dependence of /nl 

Since the sound speed can vary rapidly in our class of models, the amplitude of non-Gaussianities 
will generically display a strong scale dependence. Following [21], we define a spectral index for 
/nl as 

_ dln/NL „s 
-NG - 1 = . (8.6) 
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Figure 2: In the (csend,e) plane we plot the region corresponding to —100 < /nl < 100 with the 
choice fx = 1- The upper curve corresponds to the upper bound /nl = 100. 



In Case I, the non-Gaussian tilt arises from the ^ prefactors in (8.5). For small values of the 
sound speed, we obtain 

And since —2 < a < within the inflationary range < e < 1, we get a blue tilt for the 3-point 
function in the expanding branch, corresponding to larger non-Gaussian amplitude on smaller 



scales. As discussed in Sec. 4.2, the present range of observable scales goes from the largest 



(CMB) ones down to typical galactic sizes: i^gai/-f^cMB — 10^. We therefore expect 

/S(CMB) ^ 10=^"/S(gal). (8.8) 

Clearly we must require that non-Gaussianities are within the perturbative regime down to galactic 
scales, i.e., f^l^{gal) < 10^. In the observationally-optimistic scenario in which large-scale non- 
Gaussianities are of the order of the present experimental limits (|/nl(CMB)| ~ 100) this implies 
1*^1 ~ 2/3 (e < 1/5). Our model is less constrained if we assume instead |/nl(CMB)| = C(l), in 
which case |a| < 4/3 (e < 1/2). 

The tilt is also important in Case II. In this case, however, it cannot be approximated as 
constant over the relevant range of scales, since not all terms in (7.39) share a common {\ycnd\K)^ 



prefactor. In Fig. |3j we plot n-^c — 1 as a function of K for various parameter values. 
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Figure 3: The tilt of the bispectrum for the contracting case is plotted for different values of the 
parameters. As in the expanding case, the tilt can be large. As opposed to Case I, however, here 
the tilt can feature a remarkable running with the scale. Therefore, we plot hqtss — 1 as a function 
of K for the range of momenta between Kqmb and lO^i^cMB From top to bottom the three curves 

correspond to Cgend = 0.95, /3 = 0.1, C^end = 0.9, P = 0.01, C^end = 0.85,(3 = —0.1 



8.3 Shape of non-Gaussianities 



In light of the substantial running of the three-point function, we study the shape of the amplitude 
at fixed K, thereby disentangling shape and running effects. Following the literature [23], it is 
convenient to focus on the dimensionless ratio A/kik2k^. We will plot this quantity A/kik2k3 in 
the X2-X3 plane, where X2 = ^2/^1 and 0:3 = k^/ki. Moreover, we can restrict ourselves to the 
range 1 — X2 < < X2, where the first inequality is the triangle inequality while the second one 
avoids plotting the same configuration twice. 

Starting with the expanding case, it is instructive to first consider the limit e ^ (a — 0). In 



this regime, the last two terms in ( |7.38[ ) are subdominant and, save for the ^^3^ contribution, our 



shape reduces to the equilateral shape 

1 



A, 



equi 



OC 



i<j 



.9) 



Now, even though each term above when divided by /i;iA;2^3 is singular in the degenerate limit 
^3 — s> 0, all divergences remarkably cancel each other in the combination (8.9). One can in fact 
check that the cancellation works in the full set of amplitudes (7.38) for any value of e (a), despite 
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(c) 1^(1) - ^cqui|/(x2a;3) for a = -0.3 (d) - A,^^^/{x2X^) for a = -0.9 



Figure 4: In the top two panels we plot the non-Gaussian ampitude in the inflationary branch, 
-A^^\l,X2,x^)/{x2Xi) for Cs = 0.1, fx = 0.9, a = -0.3 (e = 0.08, top left) and a = -0.9 (e = 0.3, 
top right). Thus the dominant contribution peaks for equilateral configurations (x2 ~ 0:3 ~ 1). 
To see how this differs from DBI inflation, in the bottom two panels we subtract the dominant 
DBI-like (equilateral) contribution and plot |^*^^^ — Aeq^ui]/ {^2X3), where A^^^ and Aequi have each 
been normalized to one in the equilateral limit X2 = x^ = 1. Thus the subdominant contribution 
peaks for squashed conflgurations (x2 ~ X3 ~ 0.5). 
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other potentially divergent terms such as Xlj^j Si ^t- This is a consequence of the 

consistency relation [33] -A/ {kik2k^) — {ug — l)ki/k^ for the three-point function in the local limit 
ks <^ ki = k2, which in this case predicts vanishing amplitude since Ug = 1. Therefore, the total 
amplitude peaks in the equilateral limit ki ^ k2 ^ k^ for a generic choice of parameters, as shown 
in the top two panels of Fig. |4| and the shape is predominantly of the equilateral type. This shape 
is generally expected in models with < 1 and where higher derivatives operators are present in 
the effective action because of the shorter Jeans Length of the perturbations that easily produce 
non-linearities on shorter scales. 



equi 




0.20 



0.15 



^/-equi 0.10 



0.05 



-0.6 -0.4 -0.2 0.0 

a 



0.00 




1.0 -0.8 -0.6 -0.4 -0.2 0.0 

a 



Figure 5: The angle 9 defined in eq. (8.11 ) is plotted as a function of a in the expanding case. In 
the left panel Cg = 0.05 and fx = 0.5, in the right panel we choose instead Cg = 0.5 and fx = 0.5. 
The distinguishability between our shape and DBFs increases with |a| and with Cg. 



Nevertheless, the shape in Case I is not exactly of the form (8.9). We can unearth subdominant 



contributions by subtracting the equilateral shape (8.9) normalized at the equilateral point ki = 
/c2 = k^. The resulting amplitude, shown in the bottom two panels in Fig. |4| peaks for "squashed" 
triangles. This subdominant contribution is of course also present in standard DBI inflation, albeit 
suppressed by slow roll parameters. Because our mechanism allows for much larger values of e, 
the squashed amplitude is correspondingly more important here. 

These statements can be made more quantitative by looking at the scalar product between 
shapes Ai ■ A2, introduced by [23] 



A ■ ^2 = ^ Ai{ki, k2, k3)A2iki, k2, ks)/ (crl^al^a 



(8.10) 



hi 



The corresponding angle 612 therefore informs us on how different and distinguishable the two 
shapes are: 

^°^('^^) = (A-Ami-^2)v^- ^'-''^ 
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In Fig. [5] we plot the angle above calculated between A^^^ and ^equi as a function of the parameter 
a. As expected, the difference between the shapes increases with |a|. 




(a) -^(")(l,X2,X3)/(a;2:E3) for /3 = 0.1 (b) xa, X3)/(x2a;3) for /3 = -0.1 



Figure 6: We plot A^^^\l,X2,X3)/{x2X3) for Csend = -95, fx = ^, P = 0.1 (left) and Csend = -9, 
f3 = 0.01 (right). 



In the contracting described earlier, the curvature perturbation keeps growing outside 

the horizon, resulting in non-Gaussianities that peak in the squeezed limit, as shown in Fig. |6| 
Note, however, that the amplitude as <^ ki,k2, corresponding to squeezed configurations, is 
more singular than the local form: 



ALoc(xJ2kf. (8.12) 



This traces back to higher powers of ki in the denominator of (7.39). In particular, non-Gaussianity 
in the contracting case manifestly violates Maldacena's consistency relation, which assumes that 
C{k^) tends to a constant at long wavelengths and therefore behaves as a background wave for the 
short-wavelengths modes. This assumption does not hold in this case, while the reversed argument 
is probably the correct explanation for our shape a local type: precisely because C(^3) keeps 
growing outside the horizon the short wavelength modes effectively "feel" a different background 
evolution. The strong peak in the squeezed limit obtained here is remarkable since non-Gaussianity 
of this type usually results from entropy perturbations between two scalar degrees of freedom, such 
as in the curvaton [2] , modulated reheating [25] , pre-big bang pump field [TD] and New Ekpyrotic 
mechanisms P |29l [30l EI] . 
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9 Conclusion 



The near scale invariance of the primordial spectrum estalished by CMB and large scale structure 
observations is widely hailed as evidence for a phase of quasi de Sitter expansion in the early 
universe. In this paper we have shown that a suitably time-dependent sound speed allows for a 
much richer spectrum of early universe cosmologies, including both expanding [IT] and contracting 
backgrounds. 

In the expanding case, this mechanism was first proposed by [TTl [THl [IS], who considered a 
universe with an extremely large sound speed initially. Unlike [18] we have restricted ourselves to 
subluminal sound speed throughout. Thus the background must still be inflationary in order for 
perturbations to be stretched outside the physical horizon. Inflationary expansion of course has 
the added virtue of making the universe flat, homogeneous and isotropic. Our mechanism in this 
case alleviates the usual slow-roll conditions, from the restrictive e ^ 1/60 to e < 0.3. 

In the contracting branch, requiring that the background addresses the flatness and homo- 
geneity problems of standard big bang cosmology leads to the ekpyrotic regime e > 3 (or w > 1). 
Unlike the New Ekpyrotic mechanism which relies on two scalar fields, here a single scalar degree 
of freedom is sufficient to endow ( with a scale invariant spectrum. 

The phenomenology of non-Gaussianities is very rich and offers many distinguishing features. 
To start with, the cancelation mechanism between the sound speed and the cosmological equation 
of state does not carry over to the three-point function. Hence non-Gaussianities display a strong 
scale dependence, peaking on small scales in the inflationary branch and on large scales in the 
ekpyrotic branch. 

Moreover, the shape of the amplitude is drastically different in each case. In the expanding 
case, the amplitude is predominantly of the equilateral type, as usual with low sound speed models, 
such as DBI inflation. In the contracting case, it is predominantly of the local form, peaking for 
squeezed triangles. This contrast in shape arises because ( has different behaviors on super- 
horizon scales: whereas it tends to a constant as usual in the inflationary case, it keeps growing 
in the ekpyrotic branch. Beneath the dominant shape lies a subdominant contribution, which can 
make up as much as 10% of the total amplitude. This subdominant piece peaks for squashed and 
equilateral triangles in the inflationary and ekpyrotic branches, respectively. 

Many offspring projects naturally suggest themselves. Firstly, while our study of the non- 
Gaussianity focused on P{X, cj)) scalar lagrangians, it would be interesting to consider the more 
general effective field theory framework proposed in [56]. From this point of view, the P(X, 0) 
predictions are a subset of all possible results — hence we expect a wider spectrum of non- 
Gaussian features. Secondly, in light of the abundance of non-Gaussian features in our class 
of models, it will be enlightning to investigate the potential of near-future experiments such as 
Planck to detect some of the signatures described here, such as the strong title and the distinctive 
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shapes. Moreover, we should push the non-Gaussian exploration further and calculate higher- 
point correlation functions in this framework. The shape and tilt of the trispectrum should offer 
more ways to break degeneracies between the different cosmologies studied here. 
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A Explicit DBI Examples 

In this Section, we present explicit realizations of the scaling solutions described so far, closely 
following the analysis of [3H]. We focus on DBI lagrangians: 



Cbbi = -r^(0)v/i + /(0)(w + r'(0) - ^(0) , (A.i) 

where describes the warp factor of the throat within which the 3-brane is moving. The 
cosmological equations are 

c-^ - 1 



dlnH 1 fd(f)^^ 



dN 2M2jC, VdA 



(A.2) 



where dN = Hdt, and is given in ( 7.7[ ) 



Our goal is to determine the form of /(0) and V{(j)) corresponding to the scaling solutions of 
interest. To begin with, along the solution ~ e''^^, it is useful to think of Cs as a function of the 
classical solution for 0. Thus, 



_ dlnc,(0) d0 _ V2edc, 



where in the last step we have used (A.2). And since is assumed constant, this can be integrated 
at once for Cs(0): 

= (A.4) 

Note that we have chosen the integration constant such that the sound speed vanishes when </> = 0. 
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Similarly, on the background solution we can think of as a function of (f){N). Proceeding 
along similar lines as (A.3), e can be expressed as 

And substituting for Cs(<^) yields a differential equation for the Hubble parameter: dlnif/dln<^ = 
±2e/es, which easily integrates to 

~ 0-2^/^^ . (A.6) 

Noting that < in both Cases I and 11 — in Case 11, this assumes e > 3/2 — , we have chosen 
the sign so that \H\ decreases (increases) in the expanding (contracting) case. 

Having obtained expressions for Cs(0) and H{(f)), we can solve (A. 2) for the potential and warp 
factor: 

l/(0) = 3i^l0)M^i 

Explicitly, this gives 



= 3ff^WM^, ^1 - ^^^^^ 1 . (A.7) 



where Vq is a constant. For Case 1, where = — 2e, and in the limit of small sound speed 
(0 <^ Mpi), this gives 

^i(0) - V^o fl - II ^ . (A.9) 



3 ; 



Meanwhile, for Case 11, where = 2(3 — 2e)/5, 



2e 



10e/(2e-3) 



W)«-Vo|y-lj(^) . (A.10) 



Similarly, the warp factor can be expressed as 



which, for our explicit solutions, reduces to 

Again in the limit of small sound speed, we obtain, respectively for Case 1 and 11, 

/l(0) ~ 



Mpi 
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/ll(0) 



Ml 



PI 



(A.13) 



Thus the desired potentials and warp factors are approximately power-law in form. Note that 
in Case I we have reproduced the potential and AdS warp factor of [2S]. The exponents in V{(f)) 
and are independent of e in this case. 



B Relevant Integrals 

Below are the integrals used in the calculations of Section [7] for Case II. 
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